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Abstract 

Models with varying cosmical parameters, which were earlier regarded 
constant, are getting attention. However, different models are usually 
invoked to explain the evolution of different parameters. We argue 
that whatever physical process is responsible for the evolution of one 
parameter, should also be responsible for the evolution of others. This 
means that the different parameters are coupled together somehow. 
Based on this guiding principle, we investigate a Bianchi type I model 
with variable A and G, in which A, G and the shear parameter c 2 , all 
are coupled. It is interesting that the resulting model reduces to the 
FLRW model for large t with G approaching a constant. 
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1. Introduction 

In the past few years, evidence has mounted indicating that some constants, 
which were earlier treated as true constants, are no longer constant in cos- 
mology. The examples are — Einstein's cosmological constant A, Newton's 
gravitational constant G, the fine structure constant, etc. Different phe- 
nomenological models have been suggested to explain the evolutions of dif- 
ferent constants (let us call them parameters). However, we believe that 
there should be only one model to explain all these parameters if the under- 
lying theory is correct. Moreover, whatever physical process is responsible for 
the evolution of one parameter, should also be responsible for the evolution 
of others, implying that the different parameters are coupled together some- 
how. It should, therefore, be the evolution of the universe itself which should 
explain the dynamics of all the parameters. In this paper, we investigate 
such a model from the Einstein field equations which explains the variability 
of A, G and the anisotropy parameter a 2 simultaneously. The cosmological 
consequences of the model are also discussed. 

Now we shall describe briefly the motivation for considering the differ- 
ent parameters and their variations. The one which comes first in the list 
is undoubtedly the Einstein's cosmological parameter A, whose existence is 
favoured by the recent supernovae (SNe) la observations [1] and which is also 
consistent with the recent anisotropy measurements of the cosmic microwave 
background (CMB) made by the WMAP experiment [2]. However, there is 
a fundamental problem related with the existence of A, which has been ex- 
tensively discussed in the literature. It's value expected from the quantum 
field theory- calculations is about 120 orders of magnitude higher than that 
estimated from the observations. A phenomenological solution to this prob- 
lem is suggested by considering A as a function of time, so that it was large 
in the early universe and got reduced with the expansion of the universe [3]. 

Variation of Newton's gravitational parameter G was originally suggested 
by Dirac on the basis of his large numbers hypothesis [4]. As G couples geom- 
etry to matter, it is reasonable to consider G = Git) in an evolving universe 
when one considers A = A(t). Many extensions of general relativity with 
G = G(t) have been made ever since Dirac first considered the possibility 
of a variable G, though none of these theories has gained wide acceptance. 
However a new approach, which has been widely investigated in the past few 
years [5], is appealing. It assumes the conservation of the energy- momentum 
tensor which consequently renders G and A as coupled fields, similar to the 
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case of G in original Brans-Dicke theory. This leaves Einstein's fields equa- 
tions formally unchanged. In this context, an approach is worth mentioning 
in which the scaling of G(t) and A(t) arise from an underlying renormal- 
ization group flow near an infrared attractive fixed point [6]. The resulting 
cosmology explains the high redshift SNe la and radio sources observations 
successfully [7] . It also describes the Planck era reliably and provides a resolu- 
tion to the horizon and flatness problems of the standard cosmology without 
any unnatural fine tuning of the parameters [8]. Gravitational theories with 
variable G have also been discussed in the context of induced gravity model 
where G is generated by means of a non-vanishing vacuum expectation value 
of a scalar field [9]. Recently a constraint on the variation of G has been ob- 
tained by using WMAP and the big bang nucleosynthesis observations [10], 
which comes out as -3 x 10~ 13 yr _1 < (G/G) today < 4 x 10~ 13 yr -1 . 

Another important quantity which is supposed to be damped out in the 
course of cosmic evolution is the anisotropy of the cosmic expansion. It 
is believed that the early universe was characterized by a highly irregular 
expansion mechanism which isotropized later [11]. The level of anisotropy 
left out by the era of decoupling is only about 10~ 5 , as is revealed by the 
CMB observations. It could be that whatever mechanism diminished A to its 
present value, could have also rendered the early highly anisotropic universe 
to the present smoothed out picture. This will be our guiding principle in 
investigating the model. 

We shall keep ourselves limited to Einstein's field equations and to the 
parameters which appear explicitly therein. It would be worthwhile to men- 
tion that models with varying speed of light are recently being promoted. 
These are supported by the claims, based on the measurements of distant 
quasar absorption spectra, that the fine structure constant may have been 
smaller in the past. However, the speed of light c has a complex character 
having six different facets which come from many laws of physics that are 
a priori disconnected from the notion of light itself [12]. If it is the causal 
speed of which these theories are talking about, then one should not consider 
a varying c in general relativity unless the structure of the spacetime metric 
is changed and reinterpreted. We consider c = 1 throughout our calculations. 

We consider the Bianchi type I metric, which is the simplest anisotropic 
generalization of the flat Robertson- Walker metric and allows for different 
expansion factors in three orthogonal directions. In the comoving coordinates 
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(V = Sq), the metric can be written as 

ds 2 = -dt 2 + X 2 (t) dx 2 + Y 2 (t) dy 2 + Z 2 (t) dz 2 . (1) 

An average expansion scale factor can be defined by R{t) = (XYZ) 1 ^ 3 im- 
plying that the Hubble parameter H = R/R. 

2. Field Equations 

We consider G and A as functions of the cosmic time t. For the metric (1), 
the Einstein field equations, with perfect fluid, read 

XY YZ ZX 

XY + YZ + ZX =S ^ + A < 2 » 

X + Y + XY =-^Gwp + A (3) 
Y Z YZ 

- + - + — = -SvrG^p + A (4) 

l + f + ff = -&*Gwp + A. (5) 

Here we have assumed, as usual, an equation of state p = wp, where < 
w < 1 is a constant. The non- vanishing components of the shear tensor a^, 
defined by = u i; j + Uj-i — | gij u k . k , are obtained as 

, 4X 2 ft Z\ 



, 4 Y 2 / Z X \ 



4Z 2 [X Y\ 

- 3 z-z[x + y)- < 8) 



Thus the magnitude a 2 = a^a 10 /8 is obtained as 



x 2 r 2 z 2 /xr zx* 



(9) 
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It can be shown 2 that a 2 is proportional to R 6 , i.e., a = aR 3 , where 
a = constant. This implies that 



§--(f4+f)-^- 

Equations (2) and (9) allow to write the analogue of the Friedmann equation 

as 

3H 2 = 8nGp + a 2 + A. (11) 

So far, there has been no effect of the varying characters of G and A on the 
equations and they are formally the same as those with constant G and A. 
However, the generalized conservation equation is different from the ordinary 
one. This can be obtained either from the Bianchi identities or by using 
equations (3—5) in the differentiated form of equation (2) and can be written, 
after doing some simple algebra, in the form 

8ttG [ p + 3(l + w)Hp ] +8vrp G + A = 0. (12) 

We assume, as is common in cosmology, that the conservation of energy- 
momentum tensor of matter holds (T?j = 0) leading to 

p + 3(1 + w)Hp = 0, (13) 

leaving G and A as some kind of coupled fields: 

8vrp G + A = 0. (14) 

Equation (13 has a simple solution p = CR~ 3 ( 1+W \ where C = constant > 0. 
Equation (14) can be integrated as 



G(R) = G 



8nC 



A(R)R 3 ^ 1+w) - 3(1 +w)J A(R)R (2+3w) dR 



(15) 



where Go is a constant of integration. Equations (10—15) supply only 4 
independent equations in 5 unknowns p, R, G, A and a. In search of one 
more equation, we do some algebra in the following. 



2 By subtracting (4) from (3), and (5) from (4) and integrating the resulting equations, 
one can get — y oc XYZ , y — y oc xyz , y — ^ oc yyy ■ By squaring and adding these 



equations one gets a 2 oc 1/(XYZ) 2 
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An elimination of H between (11) and (13) gives 

t 3 (l + w) 2 ( 8 nG + - + -). 

p 3 V p p) 

Differentiating this and using (10), (13) and (14) therein, we obtain 



2^-3^ 



3(1 + w) 



1 — w 
L V 1 + w 



a 2 -A 



(16) 



(17) 



which is the central equation of our investigation whose solution will supply 
the required ansatz. Substituting (13) in (17), we obtain an equation for R 

as 



1+w 



H + 3H + 



1 — w 
1 + w 



A = 0. 



(18) 



3. Models 

If the physical processes, responsible for reducing the early highly anisotropic 
universe to a smooth present universe, are also responsible for bringing down 
the large value of A to its small present value, the two parameters a 2 and 
A must be related somehow. In view of this guiding principle, the simplest 
solution of equation (17) is 



A = 



w 



1 + w 



together with 



2^ = 3^ 



(19) 



(20) 



P P" 

Equation (19), which is our required ansatz, indicates a linear coupling be- 
tween the cosmological constant and anisotropy. The parameters G and A 
are already coupled through equation (15). We find that the model in this 
case is described by 



R = a t 2/zi - 1+w \ a = constant > 0, 
C 



P 



la 3(l+w) 

a 



t~ 



t -2/(l+w) 



(21) 
(22) 
(23) 
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G = G 



®?L „-3(l-to) 



t -2(l-w)/(l+w) _ ^25) 



AttC{\ +w) 

The model has a constant deceleration parameter q = (1 + 3u>)/2 and evolves 
to isotropy as t — > oo, with A — > and G — > G . Thus for large t, the model 
approaches the flat FLRW model which is very encouraging. It may be noted 
that though the current observations of SNe la and CMB favour accelerating 
models (q < 0), but they do not altogether rule out the decelerating ones 
which are also consistent with these observations [13]. One can even fit the 
models with zero A if one takes into account the extinction of light by the 
metallic dust ejected from the supernovae explosions [13]. 

We note that for t < 3a 2 /16nCG a 2 , G becomes negative unless w — 1 
(with Go > a 2 /8irC). One can however choose the constants a and a (which 
are arbitrary) appropriately so that G remains positive in the range of validity 
of general relativity. However, taken at the face value, the model predicts a 
repulsive gravity in the range < t < 3a 2 /lQ7rCG a 2 . For w — 1, the model 
reduces to A = 0, G = constant and a oc H = 1/31 

The model can be generalized very easily by generalizing the ansatz (19) 
in the form: 

A = 70- 2 , 7 = constant, (26) 

which can allow a negative A as well (until we have precise enough SN la 
data to rule out certain models, we should keep all our options open). Now 
with the new ansatz (26), equation (15) reduces to 

G = G - A . /T 3 ^), w^l. (27) 

47rC(l — w) 

The model starts from a big bang (or a big bang-like state) with G, | A | 
and a 2 all infinite and evolves to isotropy with G — > G and A — > as t — > oo. 
The time-evolution of R is given by 



J 3 3 V 1 — w J 



1-1/2 

dR, w^l, (28) 



where to is a constant of integration. It is hard to integrate r.h.s. of equation 
(28) for a general w 1) unless 7 = (1 — w)/(l + w) (which has already 
been investigated); or G = i.e., 7 < 0. If G = 0, equation (28) reduces to 



R 



a 



1/3 



3 1- 



w 



1 — w 



7 



1/6 



(29) 



When Go 7^ 0, A can assume both— negative as well as positive values. In 
the case of a positive A, similar kind of argument, as above, can be given 
when G becomes negative in the beginning of the universe. For a positive Go, 
equation (28) can be integrated in different phases of evolution, as is shown 
in the following. 



w 



1/3: 



1 



RV£R 2 + m - ^= sinh" 
41 




(30) 



w = : 



R 



I^ + 'o) 2 -^ 



1/3 

(31) 



where i = 8ttGG /3, m = (1 - 2 7 )a 2 /3, n = (1 - 7)a 2 /3, t = 2^/3£ and 
7 < 1/2. 

5. Conclusion 

Einstein's field equations with time- dependent G and A have been considered 
in the context of Bianchi type-I spacetime in such a way which conserves the 
energy-momentum tensor of matter. We assume that the physical processes 
responsible for the evolution of one parameter, should also be responsible 
for the evolutions of others. This means that the different parameters are 
coupled. In this view, the field equations give a trivial ansatz implying a 
coupling between a 2 (shear), A and G. The resulting model, for the baryonic 
matter, approaches the standard FLRW model in the later epochs, with G 
approaching a constant value. However, the earlier phases of the model are 
altogether different from that in the standard cosmology. For stiff matter, 
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the model reduces to the standard Bianchi type-I model with q = 2, G = 
constant, A = 0, p ~ t~ 2 and a ~ H ~ t~ x . 

Acknowledgement 

The author thanks the Abdus Salam ICTP for sending the necessary litera- 
ture whenever required under the associateship programme. 

References 

1. S. Perlmutter, et al., ApJ., 517, 565 (1999); A. G. Riess, et al, ApJ., 
560, 49 (2001); B. Narciso, et al, ApJ., 577, LI (2002); J. L. Tonry, 
ApJ. 594, 1 (2003); A. G. Riess, et al., ApJ., 607, 665 (2004). 

2. C. L. Bennett, et al., Astrophys. J. Suppl. 148, 1 (2003). 

3. Abdussattar and R. G. Vishwakarma, Pramana J. Phys. 47, 41 (1996); 
R. G. Vishwakarma, Class. Quantum Grav., 17, 3833 (2000); Gen. 
Relativ. Grav., 33, 1973 (2001); Class. Quantum Grav., 18, 1159 

(2001) ; MNRAS, 331, 776 (2002); Class. Quantum Grav., 19, 4747 

(2002) ; and the references therein. 

4. P. A. M. Dirac, Nature 139, 323 (1937). 

5. D. Kalligas, P. Wesson and C. W. F. Everitt, Gen. Rel. Grav. 24, 351 
(1992); Abdussattar and R. G. Vishwakarma, Class. Quantum Grav. 
14, 945 (1997); and the references therein. 

6. A. Bonanno and M. Reuter, Phys. Lett. B, 527, 9 (2002). 

7. E. Bentivegna, A. Bonanno and M. Reuter, JCAP, 0401, 001 (2004). 

8. A. Bonanno and M. Reuter, Phys. Rev. D, 65, 043508 (2002). 

9. A. Zee, Phys. Rev. Lett., 42, 417 (1979); L. Smolin, Nucl. Phys. B, 
160, 253 (1979); S. Adler, Phys. Rev. Lett., 44, 1567 (1980). 

10. C. J. Copi, A. N. Davis and L. M. Krauss, Phys. Rev. Lett. 92, 171301 
(2004). 

11. C. W. Misner, ApJ. 151, 431 (1968). 



9 



12. G. F. R. Ellis and J. P. Uzan, gr-qc/0305099. 

13. R. G. Vishwakarma, MNRAS, 345, 545 (2003). 



10 



